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Abstract 

The density of state for a complex N x N random matrix coupled 
to an external deterministic source is considered for a finite N, and a 
compact expression in an integral representation is obtained. 



The random matrix theory, in which the eigenvalues of random matrix 
are complex, may find some applications. For example, the two-dimensional 
electron systems under a strong magnetic field [1] or the study of neural 
network [2] is similar to the random matrix theory. 

In a long time ago, Ginibre [3] considered the complex random matrix 
theory and he obtained a density of state p{z) { z = x + iy )\ inside the circle 
in complex plane, the density of states p{z) becomes uniformly flat and is 
vanishing outside the circle. This is a generalization of Wigner's semi-circle 
law in the large N limit for the complex case. 

Recently, the random matrix theory with an external source has been 
investigated [4-16]. The external source is a deterministic, and non-random 
matrix, coupled to a random matrix. It has been discussed for a Hermitian 
random matrix [4-8] and for a chiral case [9]. Feinberg and Zee has studied 
the complex random matrix with an external source in the large N limit [10]. 
The asymmetric random matrix with external source has been also studied 
[10,11]. In the large N limit, the boundary of the density of state in the 
complex plane may be obtained by several methods. However, the expression 
for the density of state, in a finite N, is much harder for the external source 
problem. It is known that there appear interesting transitions of opening a 
gap by tuning the external source [5,13]. It may be crucial to obtain an exact 
expression for the density of state in a finite N for such problems. 

In this letter, we study a complex random matrix which couples to an 
external source matrix. We generalize the previous works for the real eigen- 
values [6,9] to this complex eigenvalue case. The density of state p{z) for the 
complex eigenvalues z is given by 



where Aj is an eigenvalue of a complex matrix M, which couples to the 
external source matrix A through the following probability distribution for 
the present case. 
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p{z) = ^ < ^tr5(a; - Ke\i)5{y - ImA^) > 
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where we write the result in which the radius of the disk is unity in the large 
N limit as a normalization. To evaluate the density of state ([^), it is useful 
to consider a chiral Hamiltonian, 



where M is a complex matrix. We denote the density of state of this chiral 
Hamiltonian by Pch{^), 

PchiX) = ^< tr5(A -H)> (5) 

where the probability distribution is P{H) = ^ exp[— A^trM'I'M]. Note that 
the eigenvalues of H are always real and appear in pairs of positive and 
negative values. Due to this chirality of the eigenvalues, the density of state 
Pc/i(A) is equal to 

p,.(A) = |A|p(A^) (6) 

where 

p(r) = ^ < tr5(r - M^M) > (7) 

in which the average distribution probability P{M) is same as Pa{M) in (0). 

As noticed by Feinberg and Zee [10], the density of state p{z) in (|l|) is 
obtained from the expression of the density of state PchW by the shift of A 
in as A A — zl . Using the well-known expressions for the complex 
delta-function 6{z) = 6{x)6{y), z = x -|- i y, 

= aA-log{zz'), (9) 

we have 

piz) = -dA' < 4trlog(^ - M){z* - Mt) > (10) 
vr iV 

Using a dispersion relation between Green function and the density of state 
Pch(A), we get 

p{z) = r dsdA4 r ^dX) 

vr Jo J-oo ts — A 

A roo i-oo \q 
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in which the external source A is shifted as A = diag{\ai — 2;|e*^% |aAr — 
2;|e*^^). In the large N limit, this p(A^) was obtained by a diagrammatic 
analysis, and the density of state p{z) was obtained by this procedure [10]. 
We consider here the finite N case, not in the large N limit, by calculating the 
chiral Pch{^) with an external source through Itzykson-Zuber integral [17]. 
A complex matrix M is decomposed as 

M = UXV (12) 

where U and V are unitary matrices and X ia a diagonal matrix. Since the 
number of real variables is 2N'^ for M, A^^ for U,V and 2N for X, we have 
2N conditions on U,V and X. It is possible to put the condition that the 
diagonal element of V is real, and X = diag{xi, ...,X]\f), Xi is real, Xi > 0. 
Note that Xi is not an eigenvalue of M, but is an eigenvalue of M'^M. Xi 
is called a singular value of M. 

Itzykson-Zuber integral for this case is known [18-20], 

dUdVe^<-^^^^^ = ^^irT}{tf}} (13) 

N\ A(a;2)A(2/2) ^ ^ 

where Y = diag{yi, ...,?/Ar) and A(x^) = YiiKji^l ~ 3;|), which is a Van der 
Monde determinant. This Itzykson-Zuber integral is obtained by applying 
a Laplacian of M to (pISl). This Laplacian reduces to the diagonal one, and 
(|13D is a zonal spherical function. 

The external source A = diag{\ai — 2;|e*^^, ...,\aN — z\e^^^) is decomposed 
as A = YU, U = diag{e^^^ , e*^^) and Y = diag{\ai — z\, {a^ — z\). This 
phase unitary matrix If can be absorbed in U. Thus, we have a diagonal 
matrix element t/i = \ai — z\ in Y . 

Using the contour representation method by Kazakov [21], and taking 
the same procedure by Brezin, Hikami and Zee (BHZ)[9], we evaluate the 
evolution operator UA{t), which is a Fourier transform of the density of state 
p(A) in d^), by noting that is an eigenvalue of M'^M, 

^^e-*^^^(t) (14) 



UA{t) = 1 < tre***^'*^ > 

N „^ N 



E f dxflx,^^det[Io{2Nx,y,)]e-''^^^+^'^- (15) 
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The coefficient of (|T^ is not important since we normalize UA{t) as Ua{0) = 
1. This expression is similar to the previous results [6,8] except that we have 
a modified Bessel function /o(2A^Xj?/j) instead of e^^'^^ as an element of the 
determinant. The modified Bessel function has an integral representation as 

/o(2v/^) = r f e^'''+--^ (16) 

Keeping the notation yf = \ai — zl"^ , we find the integral representation 
for UAit) as 

/ N 1 dd , r du 1 „-L,7)4-n 



X 



where / = ( jy-tt-jv^ )?- The contour integral over u is reduced to evaluation 
of the residue at the pole u = Nyfe~^^. The integration over q is intoduced for 
the absorption of a combinatorial factor fc!, which appears in the Xi integral 
in (0). 

One can easily find that when there is no external source ?/j = 0, the 
expression [7^=0 (^) in (O) reduces to the result of BHZ [9]. The density of 
state p(A) is given by the shift of t — > N{t + m), 

dt r du f dv , e-9-^^*^+^"^-''(^- ^ A^y? , , 



oo 



(19) 



where / = itq/{l — it), and we have made a change of variable e*^ = v/u. 
Now we consider the density of state p(z) for the complex matrix through 
(pAl) . We replace a factor as 

1 r^" 
s2 + A2 " io ' 

Inserting (|18D into ([TT]) , and doing the integration over s, A and t, we obtain 
by the change of variables, u —>■ u + a/N, q — > {l — u)q/u and a —>■ Nau and 
j3 = + 1), 

(20) 
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where the contour of the integration of v is around Nyf and the contour of 
u is around m = 1, which appears as a pole after the integration of q. 
If f{q) is a polynomial of g, we are able to prove that 

/ * r ,,^e^e-^/(,) = -e-<'-«/(.(l - (i)) (21) 
J Zm Jo u 

Thus, the integrations over q and u can be done, and we finally obtain by 
the shift V — > Nv, 

p(z) = 4-9^4/ — ^7^TT(i-^^)] (22) 

' 7tN ^Jo J 27ri /^A^ v-yf'^ ^ ^ 

where contours are taken around all yf. 

It is easy to write down the explicit form for the small N. We have 

p{z) = ie-l"^-"l' {N = 1) 

71 

1 , o 1 „-2|ai-zp _ -2|a2-zp 

= -[e-2|'^-^l' + e-2l'^-^l'-^9A*(^ ^ ^)] (iV = 2) 

TT 4 |ai — z|^ — |a2 — 

where = (oj — z)(a* — z*). It is also easy to see, when we put = 0, we 
obtain y"^ = z*z, and by the differentiation for z and z*, (B^) becomes 



Pn[z) 



If we write Ny'^ by s^, and put In{s) = Pn{z) — pm-i{z), we find In{.s) = 
/tt. This agrees with Ginibre's result (R). It is immediate to obtain 



{N-iy. 

the boundary of the density of state from (^) by the saddle point equation. 
Taking the derivative of the exponent in the large N limit by (3, and putting 
/3 = 1, which is an end point of the integral, we have as a boundary curve 
+ 1/^ = 1 for Ginibre case, and + y"^ + 2x^y^ — Sx"^ + y"^ = for the case 
when the external source eigenvalues are = ±1, N/2 times degenerated. 
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